Phase analysis of quantum oscillations in graphite 
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The quantum de Haas van Alphen (dHvA) and Shubnikov de Haas (SdH) oscillations measured 
in graphite were decomposed by pass-band filtering onto contributions from three different groups of 
carriers. Generalizing the theory of dHvA oscillations for 2D carriers with arbitrary spectrum and 
by detecting the oscillation frequencies using a method of two-dimensional phase-frequency analysis 
which we developed, we identified these carriers as (i) minority holes having a 2D parabolic massive 
spectrum p5_/2to,x, (ii) massive majority electrons with a 3D spectrum and (iii) majority holes with 
a 2D Dirac-like spectrum ±vp± which seems to be responsible for the unusual strongly-correlated 
electronic phenomena in graphite. 

PACS numbers: 81.05.Uw, 71.20.-b 



Studies of electronic properties of graphite have con- 
siderably intensified during the past decade because 
of the discovery of novel carbon-based materials such 
as fuUerenes and nanotubes constructed from wrapped 
graphite sheets 0. The quasi 2D conductivity of 
graphite occurs mostly inside the carbon layers due to 
hexagonal networks of overlapped tt— bonds. In this 
single-layer approximation the Fermi surface (FS) is re- 
duced to two points at the opposite corners of the 2D 
hexagonal Brillouin zone where the valence and conduct- 
ing bands touch each other leading to the Dirac cone 
spectrum E(j)) = ±vp±, and the charge carriers are de- 
scribed by the massless (2-1-1) dimensional Dirac fermions 

1^ Q . This point-like spectrum singularity and strong 
Coulomb coupling between fermions are assumed to be 
responsible for unusual electronic features in graphite 
such as, e. g., experimentally observed magnetic-field- 
driven metal-insulator transition 



However, to the best of our knowledge no unambigu- 
ous experimental evidence for Dirac fermions in graphite 
has been yet reported. In real graphite samples the inter- 
layer hopping leads to Pz-spectrum dispersion with open- 
ing of cigar-like FS pockets elongated along the corner 
edge H-K-H of the 3D hexaedronal Brillouin zone. The 
discussed in detail in Refs. 7^8.,M] 3D FS has a compli- 
cated multi-sheet structure and provides for the different 
groups of carriers. Band calculations show that in ad- 
dition to the two principle majority groups of electrons 
(e) and holes (h) which are located close to points K 
and H of the Brillouin zone, several minority (m) low- 
concentration groups carriers are possible. The nature 
and location of the minority pockets are very sensitive to 
the parameters of the band structure calculations and to 
the crystalline disorder. 

One can expect that the model of strongly interacting 
2D Dirac fermions is applicable to the real quasi 2D FS 
in graphite since the Dirac singularity is the topologi- 



cal property of the electronic spectra 10] that should be 
stable towards the weak 3D inter-layer coupling. 

In this Letter we make a comparative phase analy- 
sis of quantum de Haas van Alphen (dHvA) oscillations 
of the magnetization M{H) and of Shubnikov de Haas 
(SdH) oscillations of the resistance R{H) which provides 
direct evidence that the group of carriers, associated with 
majority holes (h) has the Dirac singularity in the spec- 
trum. Other groups: (e and m) have the massive spec- 
trum E{p) = p\/2m^. 

Quantum dHvA and SdH oscillations are the appro- 
priate tools to study FS properties and to distinguish 
between different types of fermion carriers. Early mea- 
surements of dHvA and SdH oscillations [tI Ii. lil ITl| 
in agreement with band structure calculations demon- 
strated that two majority (e and h) and at least one 
minority (m) group of carriers exist in graphite. 

To discriminate between normal, i. e. described by 
the massive spectrum, and Dirac fermions, we explored 
the rarely measured phase of quantum oscillations. Gen- 
erally, phase detection encounters difficulties related to 
interference of the contributions from different groups of 
carrier and to its sensitivity to errors in frequency deter- 
mination. We overcame the problem using a specially- 
developed two-dimensional phase-frequency analysis of 
its Fourier image. We distinguish also the different car- 
rier groups by applying selective pass-band filtering of 
the oscillating signal. 

The magnetoresistance R{H) and magnetization 
M{H) data were obtained on well-characterized highly 
oriented pyrolytic graphite (HOPG) sample from the 
Union Carbide Co (HOPG-UC) as described in Ref. 
5j. Briefly, low-frequency (1 Hz) and dc standard four- 
probe magnetoresistance measurements were performed 
in magnetic field — 90 kOe applied parallel to the sam- 
ple hexagonal c-axis {H || c), and at the lowest available 
temperature T — 2 K using Quantum Design PPMS-9T 
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FIG. 1: dHvA and SdH oscillations in graphite. Upper panel 
shows the region of fields 7kOe < H < 50kOe, character- 
istic for majority oscillations whereas the region of fields in 
low panel 0.25kOe < H < 2.5kOe corresponds to minority 
oscillations 

and Janis-9T magnet He-cryostats. Magnetization mea- 
surements M{H) were carried out with H \\ c by means of 
the SQUID magnetometer MPMS5 (Quantum Design). 

FigH shows the measured magnetic susceptibility x — 
dM/dH and the oscillating part of resistance Ai? (after 
substraction of the large polynomial background Ro{H)) 
as a function of the inverse magnetic field in the 
high- and low- field regions. In agreement with previous 
experiments QiSISOl' xi^^^) is a superposition of 
at least three oscillating contributions. This can be seen 
in Fig. 121 where both principal peaks mi, ei and hi and 
their second harmonic counterparts m2 62 and in spec- 
tral intensity of Fourier transformed susceptibility |x('^) 
are plotted. At the same time, only the m and e peaks 
are seen in the spectral intensity of resistance |i?(j')|. In 
other words, the SdH /i-oscillations are strongly damped. 
The corresponding oscillation frequencies i>i and their as- 
signment to the different groups of carries which we jus- 
tify below, are given in Table Note that unlike the 

widely accepted result < j/g the hole frequency 
h'h in our sample is higher then the electron frequency 

We decomposed the measured signals x(-ff~^) and 
R{H~^) onto individual m-, e- and /i-oscillations, ap- 
plying the frequency filtering with selective pass-bands 
about the 1st and 2nd harmonics of the corresponding 
resonant frequencies z/^. As shown in Fig. |3| the results 
demonstrate the generic behavior for quantum oscilla- 
tions: an initial growth of the high- field amplitude (low 
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FIG. 2: Spectral intensity of dHvA oscillations of susceptibil- 
ity |x(i')l a-nd of SdH oscillations of magnetoresistance |ii(i^)|. 
Peaks mi, 2, ei,2, hi,2 correspond to the 1st and 2nd harmonics 
of oscillations from minority electrons, majority electrons and 
majority holes. The low and high frequency plots of |-R(i^)| 
are obtained from different sets of experimental data. 



TABLE I: Frequencies i/i and phases (pi, phase factors fit, ji. 
Si and assignment of quantum oscillations in graphite. 
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Mi 


7i 




Assignment 


minority m 


3.28 





-1 


1/2 





normal h, 2D 


majority e 


46.8 


0.757r 


1 


1/2 


-1/8 


normal e, 3D 


majority h 


64.1 


TT 


-1 








Dirac h, 2D 



Landau levels) followed by a low-field Dingle attenuation 
^ ^-Ai/H ^ rpj^g low-intensity SdH h-oscillation are re- 
covered from the noisy background of R{v) around Vh- 
The sign of Ai? in Fig. O is reversed in order to re- 
cover the behavior of oscillating part of the conductivity 
Act = A(p-i) « -Ap/pl - -Ai?. 

To proceed with the phase detection, we analyze how 
the nature of the carriers influences the phase of the 
quantum oscillations, by considering the quasi 2D spec- 
trum appropriate for graphite: 

e{p) ^ e±{p±) - 2tcoap^d± pbH, (1) 

where the perpendicular dispersion s±{p±) can be either 
of the massive (parabolic) or of the Dirac (linear) type 
(see Table zt/isi? is the Zeeman splitting {fiB — 
eh/2mc). 

The ori gina l theory of dHvA oscillations of Lifshitz and 
Kosevich was developed for 3D metals with an arbi- 
trary dispersion e(p), applied to the spectrum when 
the energy spacing hiOc between Landau levels at the 
FS is smaller then the characteristic dispersion t along 
z. The other limit of almost 2D electrons was studied 
quite recently and the general expression that incorpo- 
rates both huc t and hiOc 3> t limits was derived in 
[l^j l for the case of parabolic dispersion of e± (p\). It is 
straightforward to generalize the calculations of [l3| for 
the case of arbitrary dispersion, using as in [l^ the elec- 
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FIG. 3: Quantum oscillations of susceptibility Ax{H~^) and 
of resistance AR{H~^) for different groups of carriers, ob- 
tained after two-harmonic band-pass filtering of experimental 
data. Dot-lines show the one-harmonic phase fit using the 
phase factors specified in Table H] The sign of AR is inverted 
to recover the oscillating part of the conductivity Act. 



tron orbit area S{e,pz) — np']_{e,pz) instead of ej^(pj^) 
and the Bohr-Sommerfeld semiclassical quantization: 



S{e,pz) - {n + -f)2Trh 



eH 



(2) 



where the factor < 7 < 1 is related to the topology of 
the FS: 7 = 1/2 for parabolic dispersion and 7 = in 
the Dirac case fld j. 

Reproducing the calculations of [3| in terms of 
S{e,pz) we get for the oscillating part of the magneti- 
zation: 



AM = - 



1 



S 



1 XI 



-2Trl 



(27r)2 M he dS/de ^ / sinh 



(3) 



2t 



X Jo ( 27r/- — )sin 27r/ 

nujc 



c S{e) 
eh2TrH 



7 



cos 27r/ 



TlLJc 



where T is the impurity width of the Landau level and 
parameters 

N r., .X , ""cT dS ^ eH 2ti 

S(e) = S{e,Ti2d), A = -— — « 1, = — — — 

efiH de c dS/de 

are given in Table ITll for the normal and Dirac fermions. 

Calculated with respect to the band origin (at Pz = 
0) the chemical potential e equilibrates the oscillating 



TABLE IL Spectra, Landau quantization, areas of the quasi- 
classical electronic orbits and parameters of the dHvA oscilla- 
tions for the massive (Normal) and masseless (Dirac) fermions 





Normal 


Dirac 


e(p) 


pi/2mx 


±v\p±\ 


e±{n) 


(eh/m^c) H{n+ |) 


±{2v^eh/cy/^{Hn) 


S{e) 


2nm_i_e 






eH/cm± 


ev^H/ce 


A 


2n^cmi_T/ehH 


2Tv^cTs/hev^H 



Fermi-levels of different groups of carriers and therefore 
acquires the field dependence that was shown 0| to be 
important for the very clean 2D systems with Wc 3* F and 
in the ultra-quantum limit when only low Landau levels 
{n ^ 1) are occupied. In the opposite case, we neglect 
this dependence and assume that e = ep- Note, however, 
that beats in the majority oscillations (Fig. |2Jl can be 
attributed to the conserving e " cross-talk" between e and 
h carriers We neglect also the last oscillating spin 
factor in Q since the Zeeman splitting hbH in graphite 
is much smaller than the distance between Landau levels 
hujc This tiny splitting feature is observed only in 

the high-field e-oscillations in Fig|31 

Equation includes both the 3D Lifshitz-Kosevich 
limit 12] when C = 27r/^ > 1, Jo(C) ~ 
(2/7rC)^/2cos(C - 7r/4) and pure 2D limit when C <^ 1, 
JoiC) ~ 1. In the case of 2D Dirac fermions it re- 
duces to the result obtained in Ref. [3| , whereas in 
the case of massive fermions the result of [13 is recov- 
ered. Analyzing only the oscillating part of AM, we 
find that the 1st harmonic of the magnetic susceptibility 
Ax = d{AM)/dH oscillates as: 



Axi ~ cos 



(2tt 



H 



(4) 



where the factor fi ~ sign{e) is equal to +1 for the elec- 
trons and to —1 for the holes. The topological index 7 
(see Eq. ^) is equal to 1/2 for massive fermions and 
is for Dirac fermions. The factor S reflects the curva- 
ture of the FS in the z-direction and changes from for 
a quasi-2D cylindrical FS when hojc ^ i to ±7r/8 for a 
corrugated 3D FS when huJc 3> t ( ± corresponds to the 
contribution from minimal/maximal cross section). 

In order to determine factors fj,i, 7^ and Si for each 
group of presented in Fig. oscillations we extract 
the phase of their 1st harmonics ipi, analyzing Fourier- 
transformed susceptibility x(^) i'^ the vicinity of oscil- 
lation frequencies Vi. Because the phase information 
can not be extracted from the spectral intensity plot 
|x(i^)| (Fig. EJ, we developed the method of the 2D 
phase-frequency analysis that is free from the above men- 
tioned phase/frequency uncertainty. For its illustration 
we assume that close to the resonance frequency vq the 
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FIG. 4: Contour plot of the positive part of the phase-shift 
function K{(l), v) = Ree"^x('^) for minority and majority car- 
ries. Position of maxima of K{(j), i/) determines the oscilla- 
tion frequencies Vi and phases (j>i for different groups of carri- 
ers. Upper panel presents the corresponding spectral intensity 



(i) The minority carriers are holes with a 2D massive 
parabolic spectrum. 

(ii) The majority electrons have the parabolic spec- 
trum with 3D pz -dispersion. 

(iii) The majority holes are 2D Dirac fermions. 

The identification of Dirac fermions which can be re- 
sponsible for unusual strongly-correlated electronic phe- 
nomena in graphite 0, 0| is the principal result of this 
work. 

The method proposed here for the two-dimensional 
phase-frequency analysis allows the efficient phase def- 
inition in any quantum oscillation phenomena, includ- 
ing those in low-dimensional organic conductors, in the 
mixed state of superconductors and in 2D quantum Hall 
semiconductors . 

The work was supported by FAPESP and CNPq 
Brazilian scientific agencies. I.L. thanks J.-L. Dellis for 
discussion of computation problems. 



Fourier-transformed susceptibility has the Gaussian-like 
profile x(^) ~ g'^vo^-f^ii^-i^o) ^ Constructing now the 
phase-shift function K(if, v) = Ilee~^^xi'^) that in our 
model case is equal to e"'^^''^''"' cos{(j) — 0o) we can de- 
tect both the frequency vq and the phase ipQ simultane- 
ously as the position of the maximum of K{ip, v) in the 
plane (j)-v. 

The phase-shift function for the dHvA oscillations in 
graphite is shown in Fig. 0] Determination of maxima 
of K{(p,i^) gives the collected in Table oscillation fre- 
quencies Ui and their phases ipi more precisely than the 
previous determination: (6^ = (j)f^ = 0.757r in and 
(f)h = O.TGtt, 0e = 0.647r in (l). 

By analyzing the oscillations of Act and Ax (Fig. 
and taking into account the relation Ac ^\m\H^Ax 
we conclude that the in- and out-of-phase behav- 
ior corresponds to electrons and holes with fi = +1 and 
fi — —1, respectively. For h-carriers this analysis was 
independently supported by a comparison between the 
longitudinal resistance oscillations AR and of the Hall 
resistance ARh (not shown): for holes, jminima in AR 
should correspond to minima in ARh [ijl, and this is 
what we observed. 

Knowledge of and /i^ allows the unambiguous de- 
termination of the factors 7^ (7^ = 1/2 or 0) and 6i 
{\Si\ < 1/8) that are interrelated as: 

0i = TT (sign fii - 27j/ -I- 2Si) (5) 

The analysis of the 1th harmonic parameters fii, 7^ and 
6i, given in Table |l] for each group of carriers, led us to 
the following conclusions. 
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